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ABSTRACT 

The  stress  distribution  around  a  reinforced  circular  hole  in 
a  steel  cylindrical  shell  is  determined  experimentally  for  various 
hole  sizes  or  degrees  of  reinforcement.  Results  are  compared  to 
predictions  from  flat  plate  theory,  Lure's  theory,  and  the  General 
Technology  Corporation  perturbation  technique e  Stress  concentration 
factors  based  on  the  Huber-Hencky-Von  Mises  maximum  distortion  energy 
theory  of  failure  are  computed  based  on  the  theoretical  stress  ex- 
isting in  the  undisturbed  field 0 

It  is  concluded  that  flat  plate  theory  offers  a  reasonable 
approximation  to  the  problem  of  stress  distribution  up  to  a  ratie 
of  ■   ■  ■  »  0»77>  where  a  is  the  radius  of  the  penetration,  K^   the 
mean  shell  radius,  and  t  the  thickness  of  the  shell,  all  measured 
in  inches o  The  other  theoretical  approaches  do  not  agree  well 
with  measured  values.  The  Hencky-Von  Mises  stress  at  the  periph- 
ery ef  the  hole  may  be  expected  to  become  as  great  as  three  times 
the  field  stress,  but  the  effect  of  the  hole  is  negligible  at  one 
hole  diameter  from  the  edge  of  the  hole* 
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NOTATION 

a  -  Radius  of  circular  penetration  in  cylindrical  shell. 

b0-  Radius  of  unreinforced  hole  •  Outer  radius  ef  reinforcing  plug. 

h  »  Thickness  of  reinforcing  plug, 

i  --v^T 

r  »  Radial  distance  in  (r,0)  polar  coordinate  system  centered  on  hole 
(See  figure  1,  page  ri). 

t  -  Thickness  of  cylindrical  shell  ■  R2  -  R^« 

m  -  Displacement  in  Mx"  direction 0 

r  -  Displacement  in  "y" . direction 

w  »  Displacement  of  middle  surface  of  shell 0 

(x,y)  •  Rectangular  coordinate  system  centered  on  hole  (See  figure  1,  page  ri), 

D  •  Diameter  of  circular  penetration  in  cylindrical  shell  ■  2ae 

E  -  Young's  Modulus, 

N  ■  Normal  force  per  unit  length,  subscripted  as  necessary  to  indicate 
direction,  i,e,,  N  -  force  in  "xM  direction  *C5^to 

P  •  Internal  pressure  in  cylinder,  pounds  per  square  inch, 

R,  -  Inside  radius  of  cylindrical  shell, 

R„  ■  Outside  radius  of  cylindrical  shell, 

R^  -  Mean  radius  of  cylindrical  shell  -  R^  ♦  R2 

5 

SCF  -  Stress  concentration  factor,     ^max. 

orth.0r.tieUlfield) 

where  Cy^j-,  in  the  ease  of  test  data,  is  the  maximum  measured 
stress  in  the  direction  under  consideration.  In  the  case  of  the 
theoretical  studies,  it  is  the  maximum  calculated  stress  in  the 
direction  under  consider ationc 

SS  ■  Strain  sensitivity,  microinches  per  inoh  per  psi  of  pressure. 

Vp  ■  Percentage,  by  relume,  or  reinforcement  (See  equation  /l/  on  page  32), 


<*  -  Hole  parameter  ■  — i 


«, 


V 


-ir= 


€   ■  Lineal  strain,  microi;  ches  per  inc& 

0  •  Angular  direction  in  (r,0)  polar  coordinate  system  centered  on  hole 
(See  figure  1,  page  ri), 

A  "  Azimuth  angle  of  a  meridional  plane  (See  figure  1,  page  ri). 

Z"  ■  Poisson's  Ratie. 

7r   -  3.ua$9... 

CT  ■  stress,  subscripted  as  fellows: 

Cu  -  Stress  in  nx"  direction. 

Oy    -  Stress  in  "yn  direction*, 

®e   -  Stress  in  "0"  direction , 

^r  »  Stress  in  "r"  direction, 

&i  »  <5&   ■  Principal  stresses  o 

0"HstKr  m   Hencky-Von  Mises  stress  (See  page  95  for  definition)  0 
0   ■  Displacement  f unotion 

/^  -  Complex  displacement  function 
^  •  L&Place  Operator  : 

"    "^si  +~  4r — ,    i»  rectangular  coordinates. 
D  x  *     5  y  *• 

-  j) —  +  1.    2 —  +  -I,  Js — ,  in  plane  polar  coordinates  0 
$  »  Airy  stress  function 


.  Y  - 


X 


(r,  9)  ■  Polar  coordinate  system  centered  on  the  circular  penetration 
of  the  cylindrical  shell}  r  -  radial  distance  from  the  center 
of  the  hole  measured  along  the  middle  surface  of  the  cylinder, 
6  ■  angle  between  "rn  and  a  generatrix  of  the  cylinder  passing 
through  the  center  of  the  circular  penetration. 

(x>y)  "  Rectangular  coordinate  system  centered  on  the  circular  pene- 
tration; "x"  measured  along  a  generatrix  ef  the  cylinder  pass- 
ing through  the  center  of  the  penetration,  wy"  measured  along 
a  curre  perpendicular  to  "x"  and  lying  in  the  middle  surface 
of  the  cylindrical  shell. 

SHELL  COORDINATES 


Figure  la 


SHELL  DETAILS 


Figure  lb 


COORDINATE  SYSTEMS 


Figure  1 


INTRODUCTION 

GENERAL 

The  purpose  of  this  investigation  is  to  determine  experimentally 
the  elastic  stress  distribution  around  a  radially  oriented  circular 
hole  in  a  pressurized  cylindrical  shell,  for  varying  degress  of  rein- 
forcement of  the  hole.  The  configuration  is  illustrated  in  Figures  11 
and  12,  pages  106  and  107.  This  investigation  was  suggested  by  Mr.  John 
Pulos,  Head  of  the  Fundamental  Research  Branch  of  the  Ship  Structure 
Division  of  the  David  Taylor  Model  Basin,  in  connection  with  a  program 
of  study  of  penetrations  in  cylinders  involving  multiple  holes,  different 
hole  orientations,  and  combinations  thereof. 

The  primary  goal  is  to  contrast  and  correlate  existing  theories, 
validating  them  insofar  as  is  possible  to  a  "reinforced"  case.  A  secondary 
goal  is  to  provide  empirical  data  which,  for  geometrically  similar  shapes, 
might  find  application  for  design  purposes. 
BACKGROUND 

In  recent  years,  interest  has  been  increasingly  focussed  on  the 
problem  of  relatively  large  penetrations  in  circular  shells  subject  to 
pressure.  Very  little  useful  information  exists  which  can  provide  back- 
ground and  guidance  to  the  engineer  facing  a  specific  design  problem  in 
this  field.  Although  the  problem  of  a  flat  plate  with  a  hole  is  fairly 
well  understood,  the  field  of  cylinders  with  openings  is  still  in  its 
infancy  of  investigation.  The  need  for  fundamental  research  in  this  area 
becomes  apparent  when  one  considers  submarines  going  to  ever  greater 
depths,  and  the  increasing  importance  of  nuclear  power,  with  the 
attendent  large  pressure  vessels  and  piping  penetrations  required 
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Some  theoretical  stress  analyses  of  this  oroblem  have  been 
made,  all  of  vhich  are  concerned  with  the  unreinforced  hole.  In 
addition,  experimental  verification  of  these  analyses  appears  to 
be  limited. 

Lure  (1)*,  (2)  has  given  a  solution  for  the  stress  distribution 
in  a  thin  cylindrical  shell  with  a  circular  hole.  The  results  of 
Lure's  analyses,  however,  are  valid  only  for  an  extremely  small 
hole,  one  where  ©C«  1.  The  practical  application  of  this  "pinhole" 
is  somewhat  difficult  to  comprehend. 

Savin,  in  (3),  made  use  of  Lure's  analysis  to  give  a  more 
complete  coverage  to  the  field  of  stress  concentrations  around 
holes  in  various  bodies.  Further,  in  (4),  Savin  developed  a  forata- 
lation  of  the  problem  of  a  hole  of  arbitrary  contour  in  a  thin 
shell,  but  gave  no  indication  as  to  the  solution  of  the  problem. 
It  was  stated,  however,  that  experiments  demonstrated  the  rapid 
decay  of  disturbances  in  a  uniform  stress  field  caused  by  a  circu- 
lar hole  in  the  stressed  shell,  and  that  the  effect  of  the  hole 
was  negligible  outside  of  a  distance  equal  to  the  hole  diameter, 
measured  from  the  edge  of  the  hole. 

Radok,  et  al  in  (5),  developed  an  approximate  solution  for 
the  intersection  problem,  but  as  it  is  an  energy  method  with 
approximate  boundary  conditions,  its  limitations  cannot  be  accu- 
rately predicted. 

*  Numbers  in  parentheses  refer  to  the  Bibliography  which  is  found 
on  pages  51  and  52. 
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When  a  flat  plate  with  a  small  circular  hole  is  loaded  in 
shear  by  the  action  of  two  mutually  perpendicular  forces,  equal 
in  magnitude  but  opposite  in  sign  (one  tensile,  the  other  com- 
pressive), the  shear  stress  at  the  periphery  of  the  hole  becomes 
four  times  as  great  as  the  field  shear  stress  (7),  Withum,  in  (6), 
employed  a  perturbation  technique  to  solve  the  problem  of  torsion 
of  a  cylindrical  shell  with  a  circular  hole.  Using  this  approach, 
he  showed  that  under  a  torsional  loading,  shear  stress  concentrations 
as  high  as  10  times  the  theoretical  field  shear  stress  might  be 
expected  at  the  periphery  of  the  hole.  This  is  significantly 
higher  than  the  factor  of  U   encountered  at  the  hole  in  a  flat  plate. 

Kline,  et  al,  in  (8)  using  a  similar  perturbation  method, 
described  the  state  of  stress  in  an  infinite  cylinder  under  internal 
pressure  with  a  relatively  large  circular  hole  in  the  shell. 

The  most  recent  experimental  comparison  with  theory  was  done 
by  Houghton  in  (9).  Using  a  photo-elastic  technique,  he  determined 
that  up  to  values  of  oC  -   .88,  the  stress  concentrations  given  by 
flat  plate  theory  seem  acceptable;  that  the  effect  of  shell  curvature 
was  negligible. 
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Test  Apparatus 

The  test  apparatus  consisted  of  a  pressure  vessel,  hydraulic 
pump,  flexible  tubing,  pressure  gage,  and  equipment  required  for 
the  taking  of  strain  gage  data.  The  general  arrangement  of  the 
apparatus  and  details  of  the  instrumentation  are  shown  in  figures 
2  through  4-  (pages  5  through  7).  Assembly  drawings  of  the  pressure 
vessel  and  a  detailed  description  of  the  apparatus  are  contained 
in  Appendix  E,  pages  103  through  112. 
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FIGURE  3 
Photograph  of  Assembled  Apparatus 
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DETAILS    OF   INSTRUMENTATION    OF  PRESSURE    VESSEL 

FIGURE     4 
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After  securing  the  closure  plug,  the  pressure  vessel  was  placed 
In  the  vertical  position  and  filled  vith  hydraulic  oil.  All  electrical 
leads  from  the  strain  gages  vere  attached  to  the  terminal  board  and 
the  hydraulic  pump,  piping  and  pressure  gage  were  attached  to  the 
pressure  vessel.  All  electrical  circuits  were  cheeked  for  continuity 
and  grounds. 

Pressure  was  applied  to  the  vessel  by  operating  the  hand  pump.  As 
each  desired  pressure  was  reached,  the  valve  on  the  hydraulic  line  at 
the  pressure  gage  was  secured  to  keep  the  pressure  constant  in  the 
vessel  while  strain  gage  readings  were  taken  with  the  Baldwin  Type  N 
strain  Indicator. 

All  strain  gages  were  read  at  each  pressure  by  selecting  individual 
gages  and  appropriate  dummy  gages  through  the  switching  unit  and  reading 
strains  on  the  strain  indicator.  Pressure  was  raised  in  100  PSI  increments 
to  a  maximum  of  700  PSI. 

After  all  readings  were  taken  for  a  particular  hole  opening,  the 
wires  were  disconnected,  the  hydraulic  oil  drained  out  and  the  pressure 
vessel  transported  on  a  dolly  to  the  machine  shop  where  a  larger  hole 
was  machined  out. 
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PRESENTATION  CF  RESULTS 

A.  ORIGINAL  DATA  SUMMARIES 

The  data  collected  during  the  teets  of  the  0o9$",  1.25",  l.SO",  and 

1»75"  radius  penetrations  is  presented  on  the  pages  following  in  terms  of 

t 

"strain  sensitivities'1  for  each  gage. 

These  strain  sensitivities,  in  micro-inches  per  inch  per  psi  of 
pressure,  were  obtained  by  plotting  the  individual  strain  readings  ef 
each  gage  (xcin./in.)  against  the  corresponding  pressure  (psi)  in  the 
model,  and  taking  the  slope  of  the  linear  portion  of  the  resulting  plot0 

Although  somewhat  laborious,  this  technique  has  the  advantage  of 
averaging  out  minor  errors  in  reading  the  strain  indicator  and  is  neces- 
sary to  discount  any  initial  erratic  behavior  due  to  locked-in  stresses 
in  the  shell. 

Trial  runs,  for  the  purposes  of  familiarization  with  the  apparatus 
and  checkout  of  all  associated  equipment  were  conducted  with  the  solid 
plug  in  place,  and  with  a  hole  ef  O.U75"  radius  in  the  reinforcing  plug. 
Information  received  from  the  David  Taylor  Model  Basin  indicated  that 
no  significant  increase  ef  stress  would  occur  until  at  least  one«half 
of  the  diameter  ef  the  reinforcing  plug  was  removed,  and  such  was  found 
to  be  the  case.  These  preliminary  runs  are  therefore  not  reported 0 
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TABIB  1 


SUMMARY  OF  TEST  CONDITIONS 


TEST 

DATES 

m 

a/b 

a/t 

*An 

c< 

VR(S) 

1 

15  NOV  62 

0,950 

0.500 

2.533 

0.122 

0.309 

139  oO 

2 

23  DEC  62 

1.2$0 

0.658 

3.333 

00161 

Oo537 

105  o2 

3 

22  MAR  63 

1.500 

0.790 

a.  ooo 

0.193 

0.772 

70.0 

a 

7  AFR  63 

1.750 

0.921 

U. 667 

0.225 

lo050 

28.1 
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SUMMARY  OF  TEST  DATA 
STRAIN  SENSITIVITIES  ^c  itt./\(\ ./  psc 
INSIDE  GAGES 

Gagejo*  a=0^5«  ggUaS"  ggLE  a=1.75" 

101  .260  .334  .344  .307 

102  .4.07  .472  .813  1.733 

103  .385  .460  .550  .700 

104  .385  .434  .656  1.300 

105  .222  .304  .418  .630 

106  .498  .528  .648  1.100 

107  .185  .262  .308  .500 

108  .644  .646  .662  .840 

109  .140  .140  .144  .178 

110  .626  .640  .638  .670 

111  .558  .490  .025  -.120 

112  .113  .125  .040  -1.625 

113  .684  .648  .313  .163 

114  .100  .1023  .0024  -.427 

115  .624  .617  .466  .370 

116  .127  .1074  .054  -.230 

117  .616  .652  .634  .690 

118  .122  .113  .094  -.030 

119  .580  .600  .650  .720 

120  .200  .202  .198  .200 

121  .630  .702  .240  .388 

122  .706  .842  .880  1.725 

123  .403  .445  .130  .863 

125  .634  .764  .458  .604 

124  .765  .890  .733  1.075 

126  .948  1.080  .960  1.052 

127  .584  .712  .416  .580 

128  .781  .866  .650  .863 

129  .460  .560  .278  .563 

131  .608  .690  .390  .500 

130  .783  .913  .590  .660 

132  .498  .574  .310  .450 

133  .556  .656  .294  .260 

134  .765  .844  .522  .528 

135  .461  .539  .213  .230 
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SUMMARY  OF  TEST  DATA 
STRAIN  SETTS ITIVTTIES  ^ •■  «'/».//n  / psi 


OUTSIDE  GAGES 

Gase  No. 

a=0.95w 

eU9" 

ml*l* 

KSLZS 

201 

.194 

.115 

.004 

-.500 

202 

.636 

.700 

.918 

1.338 

203 

.142 

.087. 

-  .0375 

-.470 

204. 

.618 

.660 

.778 

.980 

205 

.136 

.073 

.010 

-.300 

206 

.596 

.595 

.640 

.648 

207 

.154 

.100 

.085 

-.030 

208 

.586 

.567 

.545 

.378 

209 

.127 

.117 

.135 

.150 

210 

.586 

.583 

.546 

.430 

211 

.686 

.527 

.588 

-.600 

212 

.350 

.497 

.820 

1.420 

213 

.492 

.397 

.390 

.184 

2U 

.518 

.440 

.800 

1.438 

215 

.538 

.470 

.456 

.312 

216 

.286 

.368 

.620 

1.175 

217 

.598 

.533, 

.500 

.450 

218 

.196 

.227 

.388 

.688 

219 

.582 

.557 

.526 

.480 

220 

.160 

.150 

.168 

.230 

221 

.332 

.290 

.118 

-.355 

222 

.398 

.522 

.714 

1.300 

223 

.100 

.060 

.078 

.050 

225 

.175 

.133 

.018 

-.330 

224 

.410 

.500 

.563 

.674 

226 

.500 

.547 

.570 

.600 

227 

.175 

.143 

.044 

-.175 

228 

.444 

.490 

.476 

.488 

229 

.145 

.103 

.022 

-.200 

231 

.195 

.183 

.107 

—.020 

230 

.430 

.433 

.390 

.350 

232 

.496 

.510 

.500 

.480 

233 

.195 

.220 

.194 

.194 

234 

.450 

.473 

.425 

.350 

235 

.153 

.175 

.130 

.080 
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STFFARY  OF  TEST  DATA 
STRAIN  SENSITIVITIES  yO-in/inJp 
HOLE  INSTRUMENTATION 


S  c 


Uniaxial  Gages 
Gap;e  Location (8) 


0° 
600 
900 


Rosettes 


Data 
Considered 

Unreliable 


Gage  Looatlon(6) 


0° 

30° 
60° 
90° 


1.442 

1.140 
.437 
.115 


QftRe  Location (9) 

Leg  Angle 

agi^l 

6=1*22" 

0° 

0° 

45° 
90° 

1.740 

.478 

-.638 

3.067  « 

-.400  * 

-2.900  * 

30° 

0° 

45° 
90° 

1.490 

.507 

-.667 

3.200 

-.350 

-2.900 

60° 

OO 

45° 
90° 

.544 

-.130 
-.280 

.738 

-.763 
-.850 

90° 

0° 

45° 
90° 

.206 
.0375 
-.250 

.188 

.116 

-.688 

*  Data  considered  unreliable 
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B.  ANALYSIS  OF  DATA 

1.  PRELIMINARY 

The  strain  sensitivities  derived  for  each  gage  during  a  test 
were  reduced  to  stress  sensitivities  of  PSI  (stress)  per  PSI  (pressure). 
This  was  done  using  conventional  reduction  formulae,  sample  computations 
for  which  may  be  found  in  Appendix  D,  pages  72  through  102. 

2.  COMPARISON  OF  DATA  WITH  THEORY 

a*  Comparison  of  Lure's  and  Flat  Plate  Theory  to  Measured  Data 
(  Concentrations  of  6e   in  the  shell  ) 

A.  I.  Lure  in  (1)  described  a  method  of  computing  the 

concentrations  of  the  tangential  stress  i&e)   around  a  hole  in  a  cylinder. 

The  limiting  geometry  specified  therein  ( <X  «  1  )  was  not  strictly  met 

during  this  investigation  (  see  Table  2  ),  but  it  was  believed  that  use 

of  this  approach  for  comparative  purposes  would  not  prove  entirely  invalid. 

Accordingly,  a  set  of  stress  concentration  factors  based  on  Lure's  method 

were  computed  (  see  Appendix  B,  pages  65  through  68  )  for  selected  values 

of  0  and  r/a. 


Test 

_.  _ . —                      '  i 

Hole  Radius               oC 

1 

2 
3 

0.95"                       .309 
1.25"                       .537 
1.-50"                       .772 
1.75"                    1.050 

TABLE  2 

Using  Wang  (7)  the  solution  for  a  flat  plate  under  the  action  of 
a  biaxial  stress  field  was  found,  and  the  tangential  (6>e)  stress 
concentration  factors  were  calculated  for  selected  values  of  0  and 
r/a  (  see  Appendix  C,  pages  69  through  71  ). 
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The  tangential  stress  (&&)   sensitivities  for  the  middle  surface 
of  the  shell  were  calculated  by  averaging  the  stress  sensitivities 
of  the  inside  and  outside  tangentially-oriented  gages  at  each  gage 
location.  In  order  to  generate  the  required  stress  concentration  factor, 
each  of  these  calculated  stress  sensitivities  was  divided  by  the 
theoretical  axial  stress  sensitivity  existing  in  a  cylinder  of  finite 
wall  thickness,  as  shown  in  Appendix  D,  pages  87  and  88.  (  Although 
Novozhilov  in  (17)  considers  that  shells  whose  parameter  t/R  =  .05 
may  be  considered  as  thin  shells,  and  the  geometry  of  the  vessel  under 
test  was  such  that  t/lLj  =  .04-82,  it  was  found  that  using  the  thick- 
wall  or  Lame  stresses  resulted  in  a  closer  agreement  with  theory  ). 

The  results  of  these  computations  are  shown  on  the  following 
four  pages,  Figures  5a  through  5d,  as  dimensionless  plots  of  -jr^ 
versus  r/a  for  the  actual  and  theoretical  stress  concentrations. 
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be  Comparison  of  the  Method  of  Kline,  et  al,  to  Measured  Data 
(Concentrations  of  (^e  in  the  shell) 


As  explained  in  Appendix  A,  pages  53  through  59,  Messrs ■ 
Kline,  Dixon,  Jordan,  and  Eringen  of  the  General  Technology  Corporation, 
using  a  perturbation  technique,  described  the  theoretical  tangential  (9) 
stress  distribution  around  a  hole  in  an  infinite  cylinder  (8).  The  au- 
thors modified  this  by  the  superposition  of  the  flat  plate  solution  in 
the  longitudinal  direction  to  account  for  the  axial  stress  and  coded,  in 
FORTRAN  language,  a  program  for  the  solution  of  the  resulting  problem. 
This  was  put  inte  the  IBM  7090  computer  at  David  Taylor  Model  Basin  and 
the  results  appear  for  selected  values  of  0  and  r/a  in  Appendix  A,  pages 
63  and  6h. 

These  results  can  only  be  classed  as  disappointing.  Although  not 
entirely  unrealistic  at  the  edge  of  the  hole  (r/a  ■  1),  little  agreement 
can  be  found  between  the  actual  state  of  stress  in  the  shell  (r/a  >  1) 
and  that  predicted  by  the  General  Technology  Corporation  approach.  For 
illustrative  purposes,  two  comparative  plots  are  shown  on  figure  6  on  the 
following  pagej  one  showing  predicted  and  actual  tangential  (9)  stress 
concentration  factors  for  0  -  0°  on  the  1.5"  hole  radius  test,  and  one 
for  0  -  90*  on  the  1.75"  hole  radius  test.  Possible  reasons  for  this 
lack  of  agreement  are  advanced  in  the  section  "Discussion  of  Results" • 
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c»  Comparison  of  Theory  and  Measured  Data  at  the  Hole 
(Concentrations  ot<3&  at  r/a  -  1  for  one  quadrant; 
0  -  0°  te  90°) 


Hole  instrumentation  was  successfully  accomplished  for 
penetrations  of  1.25",  1.5%  and  1,75"  radius.  Stress  concentration  fac- 
tors were  computed  from  measured  strain  sensitivities  in  a  manner  similar 
to  that  previously  described,  the  base  being  the  theoretical  axial  stress 
in  the  cylinder. 

The  following  three  pages  compare  these  results  to  flat  plate,  Lure, 
and  General  Technology  Corporation  theory.  In  addition,  one  other  method 
is  shown;  that  used  by  Walsh  (15)  to  examine  the  stress  concentrations 
around  a  reinforced  circular  hole  in  a  submarine  hull.  Derived  from 
curves  in  Petersen  (Hi),  sample  computations  for  this  type  of  calculation 
may  be  seen  in  Appendix  D,  page  102. 
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d.  Comparison  of  Flat  Plate  Theory  to  Measured  Data 
(Concentrations  of  (5^  In  the  shell) 


The  only  available  theory  approximating  the  radial  (with 
respect  to  the  hole)  stress  distribution  in  the  shell  was  that  ef  the 
flat  plate.  Again  using  Wang  (7),  a  set  of  radial  stress  concentration 
factors  was  computed  for  selected  values  of  9  and  r/a  (see  Appendix  C, 
pages  69  through  71) 0 

In  a  manner  similar  to  that  previously  described  for  calculation  of 
tangential  (6)  stress  concentration  factors,  the  appropriate  measured 
data  was  reduced  to  radial  (r)  stress  concentration  factors,  sample 
computations  for  which  may  be  found  in  Appendix  D,  pages  92  through  9k • 

The  results  of  these  computations  are  shown  on  the  following  four 
pages,  figures  8a  through  8d,  as  dimensionless  plots  of  ~-':  versus  r/a 
for  the  actual  and  theoretical  stress  concentrations. 
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e#  Comparison  of  Measured  Data  to  the  Maximum 
Distortion  Energy  Theory  of  Failure 
(Concentration  of  the  Hencky-Von  Miaes 
Stress  in  the  Shell) 


The  Hencky-Von  Mises  stress  is  defined  (13)  ast 

where  (p'.  and  O^  are  principal  stresses.  It  is  useful  as  a  design  criterion, 
as  failure  may  be  expected  to  occur  when  the  Hencky-Von  Mises  stress  in  a 
body  reaches  the  yield  stress. 

Henoky-Von  Mises  stress  sensitivities  were  computed  for  each  gage 
location  on  the  model.  Due  to  symmetry,  the  biaxial  gages  (oriented 
along  and  perpendicular  to  the  lines  6-0°  and  0  -  90°)  could  be  as- 
sumed to  lie  in  the  directions  of  principal  stress,  while  the  data  from 
the  strain  rosettes  placed  along  the  9  «  25°  line  was  reduced  to  give 
<5,  and  C%  * 

As  is  explained  in  Appendix  D,  page  95,  the  theoretical  Hencky- 
Von  Mises  field  stress  sensitivities  for  the  inside  and  outside  surfaces 
of  the  shell  can  be  computed  from  the  Lame  (thick-wall)  stress  sensitiv- 
ities determined  previously.  This  ^y^( theoretical,  field)  then  forms 
the  basis  for  comparison,  and  the  stress  concentrations  of  the  Hencky- 
Von  Mises  stress  are  so  computed.  Sample  calculations  for  this  phase 
may  be  found  in  Appendix  D,  pages  96  through  101. 

In  order  to  show  the  effect  of  the  removal  of  reinforcing  material, 
a  constant  hole  radius  "b"  was  assumed)  that  of  the  reinforcing  plug. 
The  reinforcement  is  then  considered  to  be  distributed  evenly  around  the 
interior  ef  the  hole.  This  is,  admittedly,  somewhat  of  an  arbitrary 
assumption,  but  permits  plotting  points  on  successive  curves  one  above 
the  other  to  show  more  readily  the  effect  of  reinforcement  removal.  This 
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would  not  be  the  case  if  the  dimensionless  abcissa  were  r/a  ("a"  being 
the  radius  of  the  actual  hole)  for  then,  as  the  reinforcement  was  re- 
moved, the  hole  radius  would  also  be  enlarged,  and  the  non-dimensional 
gage  locations  would  "shift"  towards  the  hole. 

The  amount  of  reinforcing  material,  VR,  present  during  each  test  is 
defined  as  a  percentage  by  volume;  that  volume  of  reinforcement  present 
divided  by  the  amount  of  shell  material  removed  to  install  the  reinforce- 
mento  Alternately,  it  may  be  considered  as  (volume  of  the  solid  plug 
less  the  volume  of  the  hole)  divided  by  (volume  of  the  shell  material  re- 
moved to  install  the  solid  plug)c 

Computation  of  the  volume  of  the  shell  material  removed  was  not 
susceptible  to  direct  solution  by  triple  integration.  The  General  Pris- 
matoid  Theorem  from  (19)  and  Waddle' s  and  Simpson's  Rules  from  (18)  were 
used,  with  the  result  that  a  Simpson's  integration  with  half  ordinate* 
at  the  end  (the  naval  architect's  old  standby)  gave  the  best  answer;  a 
volume  of  !u2698  cubic  inches.  All  reinforcement  volumes,  Vp,  are  ex- 
pressed in  percentages  of  this  figure,  thus: 

Vfc   s    TTh   (b*  -  <Z*)  0OoJ^i£?3       Equation  /l/ 

Results  of  these  computations  are  shown  on  the  following  three  pages, 
figures  9a  through  9c,  as  dimensionless  plots  of  (dh      (measured )/ (ShY^{ the- 
oretical) versus  r/b:  each  plot  is  for  a  specific  0  value  (0  -  0°,  2J>°, 
and  90*)  and  shows  the  concentrations  of  the  Hencky-Von  Mises  stress  for 
the  inside  and  outside  surfaces  for  four  degrees  of  reinforcement. 
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f.  Variations  in  the  Direction  of  Principal  Stress 
With  Successive  Removal  of  Reinforcement 


The  direction  of  maximum  principal  strain  at  each  gage 
location  on  the  25°  leg  was  determined  by  constructing  a  Mohr's 
strain  circle  for  each  rosette.  This  direction  was,  by  definition, 
that  of  the  maximum  principal  stress.  The  directions  resulting  from 
this  procedure  are  shown  on  the  following  five  pages. 

During  the  process  of  determining  directions  of  strain,  one 
rosette  consistently  produced  a  direction  incompatible  with  the 
stress  field  of  the  shell.  As  it  was  an  inside  rosette  (  numbers 
130,  131  and  132)  immediate  direct  observation  of  the  gage  was 
not  possible,  so  all  possible  Mohr's  circles  were  drawn  to  check 
for  the  possibility  of  either  mis-orientation  of  the  gage  or  a 
mistake  in  wiring.  The  results  showed  that  the  directions  would 
be  proper  if  the  rosette  had  been  placed  so  that  gage  number  132 
were  90°  away  from  that  shown  on  the  David  Taylor  Model  Basin 
gage  placement  and  orientation  diagram  supplied  to  the  authors. 

The  end  of  the  model  with  the  wiring  feed-through  plug  was 
removed  so  that  the  gage  in  question  could  be  seen.  It  was,  in  fact, 
found  to  be  oriented  90°  away  from  the  scheduled  position.  The  gage 
placement  and  orientation  diagram  was  corrected  accordingly. 
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DISCUSSION  OF  RESULTS 

1,  Comparison  of  Measured  Data  to  Theoretical  Predictions. 

As  may  be  seen  from  Figures  5a  through  5d,  pages  16  through  19, 
the  predictions  of  tangential  (  0  )  stress  concentrations  by  Lure's 
method  are  universally  much  higher  than  those  determined  experimentally. 
For  the  cases  of  the  smaller  holes  which  more  closely  approach  Lure's 
requirement  that  °C  <<  1 ,  the  measured  data  seems  to  indicate  that 
Lure's  correction  for  curvature  in  the  body  subjected  to  stress  is 
of  doubtful  value.  Indeed,  the  theory  which  gives  the  closest  approach 
to  the  measured  data  is  that  of  the  flat  plate,  both  in  the  amount  and 
in  the  distribution  of  the  tangential  (  0  )  stress  concentrations  in 
the  shell.  Flat  plate  theory  for  these  tangential  (  9  )  stress  con- 
centrations appears  reasonably  valid  up  to  oC  =  .774-,  corresponding 
to  a  hole  radius  of  1.5".  A  similar  observation  may  be  made  for  the 
curves  drawn  at  r/a  =  1,  at  the  edge  of  the  hole  (  see  Figures  7a 
through  7c,  pages  23  through  25  ).  There,  flat  nlate  theory  and  the 
measured  data  agree  very  closely  up  to  oC  =  .774-.  (  Houghton  in  (9) 
reports  that  at  r/a  =1,  reasonable  agreement  was  maintained  up  to 
oC   =  .88  ). 

The  effect  of  curvature  upon  the  flat  plate  solution  is  evident 
when  one  examines  the  difference  between  the  0=0°  leg  and  the  9  =  90° 
leg.  At  0  =  90  ,  the  departure  from  flat  plate  theory  is  the  greatest, 
while  at  0  =  0°,  the  flat  plate  theory  predictions  are  extremely 
close  to  measured  data. 
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The  perturbation  technique  used  by  Kline,  et  al,  and  coded  by  the 
authors  for  the  IBM  7090  computer  should  not  be  entirely  discounted. 
The  oscillatory  nature  of  the  curves  (  see  Figure  6,  page  21  )  describing 
the  tangential  (  ©  )  stress  concentrations  in  the  shell  might  be 
suppressed  somewhat  by  taking  higher  order  perturbations  of  Nai  and  Nea# 
As  explained  in  Appendix  A,  page  57,  the  authors  neglected  the  second 
and  higher  order  perturbations  of  Ne/  and  Nez,  in  a  manner  similar  to 
that  done  in  (8).  One  of  the  important  steps  in  the  perturbation  technique 
employed  in  this  approach  is  the  evaluation,  at  each  non-dimensional 
radius  (  r/a  ),  of  a  series  of  Hankel  functions,  the  number  of  which 
depends  upon  the  number  of  perturbations  employed.  Since  the  Hankel 
functions  are  of  an  oscillatory  character  themselves,  it  is  suggested 
that  taking  more  terms  in  the  basic  expression  for  -g^-  would  tend  to 
bring  these  curves  more  in  line  with  the  data  and  the  other  theories. 
The  premise  that  these  additional  terms  might  make  the  theory  more 
meaningful  is  supported  by  the  behavior  of  the  predicted  tangential  (  ©  ) 
stress  concentrations  at  the  hole  where  r/a  =  1  (  see  Figures  7a  through 
7c,  pages  23  through  25  )•  Here  the  General  Technology  Corporation 
technique  behaves  in  a  more  reasonable  manner  and,  at  the  edge  of  the 
hole,  the  Hankel  functions  at  the  various  r/a  values  are  not  employed 
as  modifying  factors  of  the  stress  disturbance  set  up  by  the  hole. 

The  stress  concentrations  predicted  by  the  Petersen  (14.)  graph- 
ical  method  do  not  appear  to  be  very  meaningful.  It  seems  that  this 
approach  is  an  extrapolation  from  an  analysis  made  by  Timoshenko  of 
the  stress  distribution  around  a  beaded  hole  (  ie.,  a  small  hole  with 
local  reinforcement  in  the  form  of  a  boss  or  welded  bead  ),  and  appears 
unjustified  for  the  geometry  considered  here. 
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The  prediction  of  the  radial  (  r  )  stress  distribution  could  be 
made  only  by  flat  plate  theory.  That  this  theory  is  insufficient  to 
completely  describe  the  stresses  in  the  shell  is  evident  upon  exam- 
ination of  Figures  8a  through  8d,  pages  27  through  30.  The  effect  of 
curvature  on  the  stress  distribution  is  again  very  marked  in  the 
0  =  90°  leg.  As  in  the  case  of  the  tangential  (  9  )  stress  con- 
centrations, this  leg  shows  more  pronounced  deviation  between  actual 
and  predicted  radial  stresses  (  ^r  )•  In  addition,  the  behavior  of 
the  radial  stress  ( &  ?   )  near  the  hole  is  quite  erratic,  due  probably 
to  the  tendency  of  the  reinforcing  ring  to  distort  and  the  entire 
hole  no  longer  remaining  circular  but  becoming  elliptical  in  shape. 

2.  Comparison  of  Measured  Data  to  the  Hencky-Von  Mises  Stress 
Failure  Criterion. 

Figures  9a  through  9c  on  pages  33  through  35  indicate  that  the 
volume  of  reinforcement  has  a  marked  effect  on  the  Hencky-Von  Mises 
stress  in  the  shell.  As  the  volume  of  reinforcement,  VR,  is  pro- 
gressively reduced,  the  measured  Hencky-Von  Mises  stress  does  not 
depart  too  radially   from  that  computed  for  the  field  until  a  VR 
of  about  100$  is  reached.  Subsequently,  at  Vr  =  70$,  an  overall 
increase  of  stress  is  observed,  which  reaches  a  measured  maximum 
of  2.95  times  the  field  stress  (  inside  surface,  25°  leg  )  at 
VR  =  28.118. 
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This  finding  is  considered  quite  significant,  especially  if 
this  maximum  distortion  energy  theory  were  to  be  used  as  a  failure 
criterion.  In  a  hypothetical  case,  a  design  stress  for  a  geometrically 
similar  vessel  based  on  a  yield  factor  of  safety  of  3  might  be  con- 
sidered adequate,  and  working  pressures  assigned  accordingly.  The 
results  of  these  tests  show  that,  in  such  a  case,  the  vessel  would 
actually  be  stressed  to  a  degree  dangerously  close  to  its  yield  point. 

It  is  believed  that  between  105$  and  70$  volume  of  reinforcement, 
the  reinforcing  material  begins  to  lose  some  of  its  property  of 
rigidity  and  participates  more  completely  in  the  distortion  of  the 
hole  into  an  elliptical  hole,  rather  than  remaining  a  true  circle.  At 
a  %  of  28.1$,  the  ring  is  so  narrow  that  it  is  quite  limber  and  bends 
entirely  with  the  hole,  thus  affording  little  restraint  to  hole 
distortion. 

The  curves  also  tend  to  confirm  Savin's  (4.)  observation  that  the 
effect  of  a  hole  in  a  shell  is  negligible  at  a  distance  of  one  hole 
diameter  from  the  edge  of  the  hole.  This  distance  occurs  on  the  curves 
under  discussion  at  the  r/b  value  of  3,  and  as  may  be  seen,  the  stress 
concentration  of  the  Hencky-Yon  Mises  stress  in  this  area  are  all 
tending  to  return  to  the  theoretical  value  of  unity. 

3#  Direction  of  Principal  Stress 

The  direction  of  principal  stress  (  see  Figures  10a  through  10e 
on  pages  37  through  41  )  as  shown  by  the  strain  rosettes  at  9  =  25° 
are  of  some  interest.  They  clearly  show  that  for  small  holes  and  heavy 
reinforcement  the  direction  of  the  maximum  principal  stress  in  the 
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shell  is  not  tangent  to  the  hole  as  might  be  expected,  but  instead 
is  directed  into  the  reinforcement.  The  reinforcement  seems  to 
"soak  up"  the  normal  stresses  in  the  shell  up  to  a  VR  of  70$,  at 
which  point  the  direction  is  what  one  might  expect  to  find  in  a 
cylindrical  shell  without  a  penetration.  At  the  lower  limit  of 
reinforcement  tested,  the  direction  of  principal  stress  has  moved 
around  until  it  comes  within  8°  of  being  tangent  to  the  hole.  For 
the  hole  radii  covered  during  the  reported  tests,  the  gage  placed 
closest  to  the  reinforcing  plug  thus  has  shown  a  total  change  of 
direction  of  principal  stress  of  31 °,  due  solely  to  removal  of 
reinforcing  material. 
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CONCLUSIONS 

1.  There  is  no  theory  currently  available  which  will  completely  and 
accurately  describe  the  stress  distribution  in  a  cylindrical  shell  with 
a  circular  hole  under  the  action  of  an  internal  pressure.  The  experiments 
show  that  the  flat  plate  theory  offers  a  reasonably  close  approximation 
to  the  actual  stresses  up  to  a  value  of  ex  -  0„77»  especially  in  the  case 
of  the  tangential  (6e)  stresses.  The  radial  (or)  stresses  depart  too  rad- 
ically from  those  observed  to  permit  a  meaningful  interpretation  at  0  « 
90°,  but  are  not  unreasonable  in  the  other  areas  (9  -  0°  and  25°). 

2<>  The  behavior  ef  the  stresses  in  the  neighborhood  of  the  reinforced 
hole  is  difficult  to  evaluate .  In  the  region  of  large  reinforcement  and 
a  small  hole  (see  figure  J>a,  page  16),  the  reinforcement  tends  to  absorb 
the  shell  stresses,  and  more  than  compensate  for  the  effect  of  the  hole» 
When  the  amount  ef  reinforcement  is  decreased  below  about  100$,  the  effect 
•f  the  hole  predominates,  and  increasingly  severe  concentrations  of  stress 
result  (see  figure  5>d,  page  19  )o 

3*  The  disturbances  of  the  basic  state  of  stress  in  the  shell  decay 
rapidly  with  distance  from  the  hole  and  at  a  distance  of  one  hole  diam- 
eter from  the  edge  of  the  hole,  the  effect  of  the  hole  seems  negligible. 

U.  Assuming  the  maximum  distortion  energy  theory  of  failure,  the  most 
critical  area  appears  to  be  the  inside  surface  of  the  shell  around  9  - 
25°.  This  seems  of  interest,  because  the  strain  rosettes  were  placed 
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along  this  line  as  a  result  of  photeelastie  studios  mado  on  the  sano 
model  at  the  Darld  Taylor  Model  Baein,  using  the  Zandtmann  birefringent 
eeatlng  technique.  The  photoelastie  study  was  necessarily  condueted 
only  on  the  outside  surface,  and  the  measured  strain  data  indie ates  that 
on  the  outside  surface,  the  most  stressed  area  is  that  at  0  -  90° c  Am 
explanation  of  this  might  be  seen  when  the  geometric  considerations  are 
examined.  The  reinforcing  plmg  thickness  was  such  that  at  8  •  0°,  it 
was  tangent  to  the  outer  surface  of  the  shell,  and  at  0  -  90°,  it  was 
tangent  to  the  inner  surfaee.  Therefore  relatirely  sharp  structural 
discontinuities  exist  en  the  outside  surface  at  0  ■  90°,  and  on  the  in- 
side surface  at  8  -  0°o  The  measured  maximum  Hencky-Vcn  Mises  stress 
concentration  is,  in  fact,  found  on  the  outside  surface  at  0  -  90°,  yet 
is  seen  at  0  -  25°  on  the  inside.  The  difference,  however,  between  that 
at  0  -  0°  and  8  -  2$°   is  less  than  2%.     Therefore,  in  addition  to  the 
effeet  of  this  hole  in  producing  a  stress  concentration,  the  effect  of 
the  plug  itself  in  producing  discontinuities  in  the  structure  may  be 
considerable o 

It  can  be  concluded  that,  with  the  vessel  geometry  under  consider- 
ation, a  concentration  ef  almost  three  times  the  Hencky-Von  Mises  theoret 

leal  field  stress  may  be  expected  to  occur  in  similar  circumstances 
for  the  inside  surface  at  0  -  0°  and  25°  (see  figures  9a  and  9b,  pages 
33  and  3U). 

« 

£.  Despite  the  fact  that  the  cylinder  was  stress-relieved  both  after 
initial  fabrication  and  after  installation  of  the  reinforcing  plug, 
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loeked-in  stresses  were  believed  to  exist  in  the  shell*  This  obser- 
ration  is  made  based  on  the  behavior  of  the  strain  sensitivities  as  read 
from  the  pressure-strain  plots*  The  first  two  of  the  reported  tests  were 
based  on  data  taken  during  runs  conducted  between  0  and  500  psi  internal 
pressure.  The  linear  portions  of  the  resulting  plots  almost  always 
failed  to  include  the  zero  strain,  zero  pressure  point.  The  last  two 
reported  tests  used  a  zero  strain  point  arbitrarily  taken  at  200  psi 
initial  internal  pressure,  and  the  pressure-strain  plots  were  made  up 
to  and  including  700  psi.  The  plots  that  resulted  from  this  technique 
were  much  more  "well  behaved",  and  tended  to  give  one  more  confidence  in 
their  use„  This  technique  of  disregarding  initial  readings  and  conduct- 
ing tests  at  somewhat  higher  leads  is  commended  to  future  investigators • 
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RECOMMENDATIONS 

1,  A  test  series  on  a  similar  model,  utilizing  the  same  hole 
sizes  as  studied  herein,  but  with  a  thicker  reinforcing 
plug  could  be  made.  This  series  would  be  able  to  provide 
more  information  on  the  effect  of  volume  of  reinforcement, 
V«,  on  the  Hencky-Von  Mises  stress  in  the  shell. 

2.  Set  up  the  computations  and  programming  necessary  for  second 
and  higher  orders  of  perturbations  used  in  the  General 
Technology  Corporation  scheme.  It  is  believed  that  this 
will  cause  the  stress  predictions  of  this  method  to  agree 
more  closely  with  those  actually  measured  in  tests. 

3»  When  utilizing  electrical  resistance  strain  gages,  always 
use  data  derived  from  the  linear  portion  of  a  load  versus 
strain  plot.  In  this  way,  erratic  gage  behavior  due  to  any 
initial  stress  in  the  model  does  not  affect  the  results. 
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APPT7TDTX  A 

Derivation  of  a  PORTRAIT  Program  for  the  IB!'  7090  Computer 
Usinr:  tha  General  Technology  Corporation  1  Method  of  Cpmnutinp 
Stress  Concentrations  as  lodified  by  the  Authors. 

The  method  used  in  this  analysis  represents  the  approach  of  T.escrs* 
L#V.  Kline,  H.C.  Dixon,  H.F.  Jordan  and  A.C.  Erinfen,  whose  Technical 
Report  Number  3-1  of  the  General  Technology  Corporation,  entitled 
"Stresses  in  Pressurized  Cylindrical  ^hell  >.Tith  Circular  Cutout",  has 
been  previously  cited  as  reference  (8).  The  follox/ing  is  meant  to  act 
as  a  brief  description  of  their  method,  and  to  show  how  a  flat  plate 
solution  for"  the  uniaxial  stress  condition  was  superimposed  ur>on  the 
problem  to  account  for  a  finite  cylinder. 

A  simplification  of  the  problem  description  was  initially  achieved 
by  considering  a  flat  strip  containing  a  circular  hole  which  is  used 
to  construct  a  cylindrical  shell  as  is  shown  below: 


*X 


(a) 


(b) 
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Although  the  hole  geometry  is  not  precisely  the  same  as  that 
obtained  by  piercing  the  shell  with  a  right  circular  cylinder  of 
radius  "a"  whose  axis  is  normal  to  the  axis  of  the  main  cylinder, 
the  difference  between  the  two  closed  space  curves  is  small  for 
relatively  small  hole/cylinder  ratios  (^  ^  ~+  )  .  It  should  be  noted 
that  the  coordinate  lines  "T"  and  "r"  remain  in  the  curved  surface. 

If  no  hole  is  present  in  the  shell,  the  hoop  stress,  £y  ,  due 

to  an  internal  pressure  "P"  is  given  by  the  "pressure  vessel"  formulat 

/  J   =  PR 

The  radial  stress,  £r,  and  the  tangential  stress,  ^©,  at  any 
point  on  the  circle  of  radius  "a"  can  then  be  described  by  considering 
the  conditions  for  equilibrium  of  an  extremely  small  triangular  element 
on  the  edge  of  the  circlet 


£r=  fy  ( \  -  cos  ze>)  =  ££(>-cosze) 

6e  =  4*    sin  IB   =    ||-    rin  2  6 

By  cutting  a  circular  hole  in  the  shell  and  applying  the  negative 
of  the  above  stresses,  equilibrium  is  maintained.  Therefore,  the  prob- 
lem of  determining  the  stresses  in  a  pressurized  shell  with  a  circular 
cutout  was  reduced  by  (8)  to; 
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(a)  Determining  the  stresses  in  the  unpressurized  shell  due 
to  the  effects  of  -  ^.and  -^acting  on  the  edge  of  the 
hole,  and 

(b)  Superposing  these  stresses  determined  in  (a)  upon  those 
due  to  pressure  in  the  cylindrical  shell  without  a  hole. 
Upon  superposition,  the  tiro  sets  of  stresses  mist  then 
elirinnte  one  another  at  the  hole. 

The  main  problem,  that  of  (a)  above,  required  that  the  ar^ronriate 
differential  equations  of  cylindrical  shells,  expressed  in  cylindrical 
coordinates  (r,©),  be  solved.  To  this  end,  use  vas  made  of  the  complex 
displacement  function- 

f  =  w  h-  L  ")L  <£ 

where,  w  =  displacement  of  the  middle  surface  of  the  shell 

i  *vcr" 

£  =  Airy  Stress  Function 

X-efc£(i-f*JJ* 

E  -  Young's  I'odulus 
~V~  ~   Foisson's  Ratio 
Expressed  in  cylindrical  coordinates  (r,©),  the  differential  equation 
for  V  is: 


A' 


3(/~  3^Z) 
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The  solutions  of  the  equations  were  subject  to  boundary  conditions 
such  that  at  the  edge  of  the  holo  the  stresses  be  -  ^rand  -  6&9   while 
bending  moment  and  normal  shear  be  zero,  and  that  the  stresses  vanish 
as  nr"  approaches  ®o# 

The  quantity  to  be  determined  in  the  approach  was  the  normal  force 
in  the  9  direction,  or  He,  which  can  be  characterized  as: 


ll©  25   _ _— -~^- 


s 


where,  £  =  fl  r 


To  effect  this,  the  unknown  function  was  expanded  into  a  Fourier 

series,  asi 

He  (r,6)  =  ^  M«n(^)co5  2n6 
r\  ~  O 

This  reduced  the  problem  by  the  separation  of  variables  to  the 

solution  of  differential  equations  for  the  Fourier  expansion  functions 

Hen  (r),  but  the  equations  still  were  not  able  to  be  directly  solved. 

A  perturbation  scheme,  similar  to  that  employed  by  T.Tithum  in  (  6  ),  was 

next  used,  wherein  the  Fourier  functions  were  expanded  as: 

Hen  (r)  «^  ^jl(r)   =  i\ll°<t(r)   4-  hi&\  (r)  +  K|£*(r)  +  ' 

<°>,  i 

The  zero  order  perturbation,  Hen  (r),  T./as  the  solution  for  the 

limiting  case  of  a  flat  plate,  while  the  higher  order  perturbations 
gave  the  appropriate  corrections  for  the  curvature  of  the  cylindrical 
shell. 

Solutions  of  the  corresponding  differential  equations  for  each 
Nen  (r)  wore  found  as  polynorials  containing  a  bounded  number  of 
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undetermined  constant  coefficients.  Recursion  relationships  were 
established  which,  with  the  Ixrandary  conditions  previously  described, 
determined  the  forces  Hen  (r)  completely. 

It  was  felt  that  the  second  and  higher  order  perturbations  of 
He1  were  quite  small  and  could  be  disregarded,  so  that  the  resulting 
force  per  unit  length  due  to  the  stress  loading  at  the  hole  was: 

"«=  N^t-lN'eT  +  Nei'j  cos  2  3 

In  order  to  obtain  the  solution  of  the  problem  of  a  pressurised 
cylinder  with  a  hole,  the  force  due  to  pressure  in  the  cylinder  without 
a  cutout  was  added,  which  can  be  expressed  as  -r^-  (  |  -+•  cd-f  26  j 
where  Ny   =  6y  t.  Therefore, 

n!=  m0  +  ^.(i+  cos  ze) 

to)    ,    a/.\  ,    f  ilic)  .    ./CO  ^    N* 


—    t 


(Hi°l  *J£)*-(u.l)+Ul*:  +  }&)  -os  ze 

This  completed  the  study  of  Kline,  et  al,  but,  as  may  be  seen,  it 
neglected  the  effect  of  any  andal  stress,  or^x,  in  the  direction  of 
©  =  0°.  To  eorreet  for  this,  the  case  of  a  flat  plate  with  a  hole  under 
uniaxial  stress  of  the  proper  magnitude  (^x.  =  -£-)   was  superposed  upon 
the  Ne  previously  determined,  by  the  follox/ing  procedure: 

Expressing  the  flat  plate  solution  (  7 )  in  the  notation  employed  by 
Kline,  et  alj 

where,  Nx  =  ^x  t 
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This  is  added  to  He  to  givej 
ITeT  ss  HeTotal  =  He  +  II© 

It  was  desired  to  transform  the  equation  into  one  yielding  a 
non-dimensional  stress  concentration  factor,  based  upon  the  nominal 
stress  e:dsting  in  the  axial  direction  (  the  lowest  stress  ).  Dividing 
both  sides  of  the  TTeT  equation  lyyNy  *a6yt$   and  realizing  that  6  %   =  ~*-t 
wo  have  an  expression  for  the  stress  concentration  factor,  or 

.     (o\  j_  \  i     (°)    ft)       <?       \  ^     y 

Stress  Cone.  Factor  =r£^~rU  F1  •  'V"*"  ^^T^^t  ^T^'^bJC03  {r~ 

Sxpressions  for  the  Hen'  terms  were  derived  from  reference  (g  ), 
and  one  case,  that  of  the  smallest  penetration,  was  worked  out  by  hand 
with  the  aid  of  reference  (t))#  A  program  for  the  IBM  7090  coin-outer 
was  subsequently  compiled,  and  the  hand  calculations  used  in  "de-bugging* 
the  initial  runs*  A  few  rl stakes  were  found  and  corrected,  and  the  whole 
program  re-run,  varying  the  hole  sizes  in  the  input  information  to  cover 
those  cases  vhich  were  experimentally  investigated.  A  copy  of  the  coT-nuter 
program  in  FORTRAN  language  and  the  results  derived  from  it  appear  on  the 
following  pages. 
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The  following  restrictions  apply  vhen  using  the  General 
Technology  Corporation  approach  I 

(a)  If  R*)£<C/£f  the  shell  is  too  thick  for  application 

of  thin  shell  theory. 
W  If  ViR'i'?' *  </*  4>  tile  Perturbation  technique  does  not 

produce  rapid  convergence,  and  more  terms  than 

Neo   ,  hl(£?    and    Me/       ***  necessary  to  eenpute 

the  stress  concentrations. 
(o)  If   /Rt^yz  §  the  circle  and  the  space  curve  obtained 

by  intersecting  two  cylinders  begin  to  differ  appreciably. 
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FOR  XRHP  10 
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PRINT2  XRHP   50 

SQ2=SQRTF(2. I  XRHP   60 
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LO=XLOCF(BES( 1 > I  XRHP  170 
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ZC1K  I  >  =  BPH«X(  1  I/SQ2  XRHP  280 

ZB12I I  >  =  2.*ZB11(  I  l+X(  I  )»«2»ZC10(  I  I  XRHP  290 

ZC12I I  >=2.«ZC11(  I l-XI  I >»»2*ZB10<  I  1  XRHP  300 

ZB13I I >=X( I )«»3/SQ2»BPH  XRHP  310 

ZC13I1  )=-XI  1  >«3/SQ2»BMH  XRHP  320 

All. 11=0.  XRHP  330 

All. 21=1. /Xll)*»2  XRHP  340 

Al 1.3I=ZC10I1 I  XRHP  350 

A(  1.4) =Z610I 1  I  XRHP  360 

AI2.1I-0.  XRHP  370 

AI2.2)=-2.*AI 1,21  XRHP  380 

AI2.3)=-2.»AI1.3)+ZC11I 1)  XRHP  390 

AI2.4I=-2.«AI 1.4I+Z811I1)  XRHP  400 

AI3.1>=4.2«AI 1.2)  XRHP  410 

AI3.2I-0.  XRHP  420 

AI3. 31=4. 2»AI 1.41-2. 7»ZB11 (1  I+ZB12I  1  I  XRHP  430 

AI3.4I=-4.2»AI1.3I+2.7*ZC11I 1 I-ZC12I 1)  XRHP  440 

AI4,1)=A(3.1]  XRHP  450 

AI4.2)=0.  XRHP  460 

A(4.3)=4.2»AI  1.41+3. 9»ZBlllll-3. 'ZB12I1  I+ZB13I1)  XRHP  470 

AI4.4)=-4.2«AI 1.31-3. 9«ZC11I 1 1+3. «ZC12I1 I-ZC13I 1)  XRHP  480 

ERASE.R  XRHP  490 

R(l.l)=-Xlll**2/4.  XRHP  500 

R(1.2)=R(1.1I  XRHP  510 

R(2.1)=2.«R( 1.1)  .^  XRHP  520 

0ETR=1.  ~b"~  XRHP0525 

XX=XSIME0F(4.4.2.A,R,DETR.BES)  XRHP  530 


JOB  059217XR 

LOAD 

TRA 


XRHP  PR1-9S8 


LOAOING  LIB 

TRANSFER  VECTORS 
0SSH1  00000 
01611  00002 


05160  00000 
OLSMU  00000 


05337  00000 
01523  00000 


SUBROUTINE  ENTRY  LOCATIONS 


05167  00000 
03651  00000 


ENOJOB 

05SU1 

SYSTEM 

05512 

CLKOUT 

05513 

XLOC 

05160 

L0G10 

05202 

LOGE 

05211 

LONE 

05167 

SORT 

05105 

COS 

01661 

COSO 

01617 

EXPA 

01616 

EXPI3 

01616 

(RLR) 

01536 

(TSBI 

01537 

(MLR) 

01510 

(FPTI 

01511 

(SPHI 

01523 

(STH) 

01521 

(CSH) 

01525 

(TSHI 

01526 

(EFT) 

01531 

(RWTI 

01532 

(FID 

01533 

(RTN) 

01531 

6ESJ 

01701 

BESY 

01701 

BESI 

01707 

8ESK 

01712 

05105  00000 
01701  00000 


EXP  05337 
SIN  01666 
EXP(2  01511 
XXXRET  01512 
(SCHI  01527 
XSIHEQ  01275 
000000  00117 


01617 

00000 

00100 

00017 

LOG 

05213 

SIND 

01651 

(STBI 

0153  5 

(ULTI 

01513 

(BSTI 

01530 

XOETRM  01271 


STRESS  CONCENTRATION  FACTORS 
METHOD  OF   GEN.  TECH.  ♦  F.  P. 


HAWK  ANO  PIERCE 


CASE  NO 

1-1 

CYLI 

JDER  PARAME 

TERS 

R  =   7.780000     T" 

0.375000 

NU=0.300 

A«   0 

475000 

RADIUS 

RANGE 

1.0000 

0. 
7.1814 

10.00 
6.9292 

20.00 
6.2031 

THETA  RANGE 
30.00      45.00 
5.0907     3.0000 

60.00 
0.9093 

75.00 
-0.6212 

90.00 
-1. 1814 

1 . 1000 

7.5450 

7.2552 

6.4208 

5.1423 

2.7397 

0.3370 

-1 .4219 

-2.0657 

1 .2500 

7.3590 

7.0635 

6.2128 

4.9095 

2.4600 

0.0105 

-1.7826 

-2.4390 

1 .5000 

6.5186 

6.2562 

5.5005 

4.3427 

2.1667 

-0.0093 

-1.6023 

-2. 1853 

1 .7500 

5.6894 

5.4663 

4.8239 

3.8396 

1.9898 

0.1400 

-1.2142 

-  1 . 7098 

2.0000 

5.0 188 

4.8292 

4.2833 

3.4469 

1.8750 

0.30  31 

-0.8476 

-1 .2688 

3.0000 

3.5179 

3.4063 

3.0848 

2.5923 

1.6667 

0.7410 

0.0634 

-0. 1846 

4.0000 

2.8817 

2.8041 

2.5804 

2.2377 

1.5937 

0.9498 

0.4783 

0.3058 

5.0000 

2. 541 1 

2.4819 

2.3116 

2.0506 

1.5600 

1.0694 

0.7103 

0.5789 

10.00  00 

1.7880 

1.7715 

1.7241 

1.6515 

1.5150 

1 .3785 

1.2786 

1.2420 

CASE  NO 

.   1-2 

CVL1N0ER  PARAMETERS 

R=   7.780000     T= 

0.375000 

NU=0.300 

A=   0 

.950000 

RADIUS 
RANGE 

0. 

10.00 

20.00 

THETA  RANGE 
30.00      4  5.00 

60.00 

75.00 

90.00 

1 .0000 

5.3222 

5. 1821 

4.7789 

4.1611 

3.0000 

1.8389 

0.9889 

0.6778 

1 . 1000 

5.7161 

5.5366 

5.0198 

4.2279 

2.7397 

1.2514 

0. 1620 

-0.2368 

1 .2500 

5.69  14 

5.4965 

4.9354 

4.0757 

2.4600 

0.8443 

-0.3385 

-0.7714 

1 .5000 

5.1553 

4.9750 

4.4561 

3.6610 

2. 1667 

0.6724 

-0.4216 

-0.8220 

1 . 7500 

4.5647 

4.4094 

3.9623 

3.2772 

1.9898 

0.7024 

-0.2401 

-0.5851 

2.0000 

4.0585 

3.9268 

3.5477 

2.9668 

1.8750 

0.7832 

-0.0160 

-0. 3085 

3.0000 

2.7418 

2.6770 

2.4903 

2.2042 

1.6667 

1.1291 

0.7356 

0 .  59 1 5 

4.0000 

1.9307 

1.9104 

1.8519 

1.7622 

1.5937 

1.4253 

1.3019 

1 . 2568 

5.0000 

1.3506 

1.3632 

1.3996 

1.4553 

1.5600 

1.6647 

1 .7413 

1 . 7694 

10.0000 

2.3461 

2.2960 

2.1516 

1.9305 

1.5150 

1.0995 

0.7953 

0.6839 

CASE  NC 

.   1-3 

CYLINOER  PARAMETERS 

R=   7.780000     T= 

0.375000 

NU=0.300     A= 

.250000 

RADIUS 
RANGE 

0. 

10.00 

20.00 

THETA  RANGE 
30.00      45.00 

60.00 

75.00 

90.00 

1 .0000 

4.7766 

4.6694 

4.3609 

3.8863 

3.0000 

2.1117 

1.4615 

1  .2234 

1 . 1000 

5.0569 

4.9171 

4.5148 

3.8983 

2.7397 

1.5811 

0.7329 

0.4225 

1  .2500 

4.9807 

4.8287 

4.3910 

3.7203 

2.4600 

1 .1997 

0.2770 

-0.0607 

1  .5000 

4.4541 

4.3162 

3.9189 

3.3104 

2. 1667 

1.0229 

0.  1857 

-0. 1208 

1 .7500 

3.8790 

3.7651 

3.4370 

2.9344 

1.9898 

1.0452 

0.3537 

0. 1006 

2.0000 

3.3641 

3.2743 

3.0157 

2.6195 

1.8750 

1.1305 

0.5854 

0. 3859 

3.0000 

1.8192 

1.8100 

1.7835 

1.7429 

1.6667 

1.5904 

1.5346 

1.5141 

4.0000 

0.8334 

0.8792 

1.0113 

1.2136 

1.5937 

1.9739 

2.2523 

2. 354 1 

5.0000 

0.5463 

0.6074 

0.7834 

1.0531 

1.5600 

2.0669 

2.4379 

2.5737 

10.0000 

2.4432 

2.3872 

2.2261 

1.9791 

1.5150 

1 .0509 

0.7111 

0. 5868 

CASE 

NO 

.   1-4 

CYLINDER  PARAMETERS 

R  =   7.780000     T= 

0.375000 

NU=0.300 

A  =   1 

.500000 

RADIUS 

THETA  RANGE 

RANGE 

0. 

10.00 

20.00 

30.00 

45.00 

60.00 

75.00 

90.00 

1.0000 

4.6551 

4.5553 

4.2679 

3.8275 

3.0000 

2. 1725 

1.5667 

1.3449 

1.1000 

U. 7769 

4.6540 

4.3003 

3.7583 

2.7397 

1.7211 

0.9754 

0.7025 

1.2500 

"1.5737 

4.4462 

4.0792 

3.5166 

2.4600 

1.4032 

0.6295 

0.3463 

1.5000 

3.9387 

3.8319 

3.5242 

3.0527 

2.1667 

1.2606 

0.6320 

0.3946 

1.7500 

3.2800 

3.2022 

2.9782 

2.6349 

1.9898 

1.3447 

0.8724 

0.6996 

2.0000 

2.6761 

2.6278 

2.4887 

2.2756 

1.8750 

1.4744 

1. 1812 

1.0739 

3.0000 

0.8187 

0.8699 

1.0171 

1.2427 

1.6667 

2.0906 

2.4010 

2.5146 

<t. 0000 

0.09411 

0.1846 

0.4452 

0.8441 

1.5937 

2.3434 

2.8922 

3.0931 

5.0000 

0.8219 

0.6664 

0.9945 

1.1909 

1.5600 

1.9291 

2.1993 

2.2961 

10.0000 

0.9742 

1.0068 

1.1007 

1.2446 

1.5150 

1.7854 

1.9834 

2.0558 

CASE 

NO 

.   1-5 

CYLINDER  PARAMETERS 

R=   7.780000    T  = 

0.375000 

NU=0.300 

A  = 

.750000 

RADIUS 

THETA  RANGE 

RANGE 

0. 

10.00 

20.00 

30.00 

45.00 

60.00 

75.00 

90.00 

1.0000 

4.7749 

4.6679 

4.3597 

3.8875 

3.0000 

2.1125 

1.4629 

1.2251 

1.1000 

4.6878 

4.5703 

4.2321 

3.7138 

2.7397 

1.7656 

1.0525 

0.7915 

1.2500 

4.2856 

4.1755 

3.8585 

3.3728 

2.4600 

1.5472 

0.8790 

0.6344 

1.5000 

3.4352 

3.3567 

3. 1384 

2.8009 

2.1667 

1.5324 

1.0681 

0.8982 

1.7600 

2.6026 

2.5657 

2.4593 

2.2962 

1.9898 

1.6634 

1.4591 

1.3770 

2.0000 

1.8411 

1.8431 

1.8490 

1.8580 

1.8750 

1.8920 

1.9044 

1.9089 

3.0000 

-0.2320 

-0.1175 

0.2122 

0.7173 

1.6667 

2.6160 

3.3110 

3.5654 

14.0000 

0.0275 

0. 1220 

0.3940 

0.8106 

1.5937 

2.3769 

2.9501 

3. 1600 

5.0000 

2.1596 

2. 1234 

2.0193 

1.8598 

1.5600 

1.2602 

1.0407 

0.9604 

10.0000 

2.1393 

2.1017 

1.9933 

1.8272 

1.5150 

1.2028 

0.9743 

0.8907 

APPENDIX  B 

tolTfffff  fi£  ifaft  Stress  Concentration  ground 
£  Hole  by,  the.  Method  fi£  A*.  Lt  fettj 

A  method  of  anrlytieally  determining  the  stress  concentrations 
In  the  neighborhood  of  a  hole  In  a  circular  cylinder  was  formulated 
by  A*  I.  Lure',  the  derivation  of  vhioh  appears  in  ( 1  )$   ( 2 )  and 
(  3  )•  Of  these,  the  last  mentioned  is  particularly  poor;  numerous 
omissions,  erroneous  inclusions,  and  undefined  assumptions  were 
noted*  This,  however,  was  felt  to  be  the  fault  of  the  translator, 
not  of  the  derivation  proper* 

The  equations  initially  used  by  Lure'  were  approximate  in  nature, 
taking  into  effect  only  the  *w*  or  radial  component  of  displacement 
and  neglecting  the  "u"  and  "v*  components  when  calculating  changes 
in  curvature  and  twist*  This  approximation  reduced  the  problem  to 
the  determination  of  a  displacement  function  <j>  vhioh,  In  the  absence 
of  surface  loads  on  the  shell,  satisfied  the  differential  equation} 

where,  A  *  -ttt  -t*  t*. 

x  s  abclssa  along  a  generatrix  of  the  cylinder 

y  st  k  X.9  Xbeing  the  asrfjmith  angle  of  a  meridional  plane 
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By  using  a  complex  stress  function,  expressions  for  the  forces 
and  moments  in  the  cylindrical  shell  were  determined,  which  contained 
a  number  of  undetermined  constants.  At  this  moint,  a  further  restriction 
was  placed  on  the  solution j 


c<   J-± 


TZ^K  V  r 


M 


This  assumption  was  made  "r.  order  to  be  able  to  neglect  terms  on 
the  order  of  £  tt   "^^  restricts  the  application  of  the  solution, 
essentially,  to  a  pinhole  in  a  thin  shell. 

Lure  nerrb  applied  a  cylindrical  coordinate  system  to  the  shell, 
centered  about  the  hole,  such  that  :•:  =  r  cos  *  and  y  =  r  sin  €",  whore 
the  curves  9  =  constant  would  be  helices  when  in  the  shell  or  radial 
lines  terminating  at  the  origin  if  the  cylinder  were  developed  into 
a  plane. 

The  stresses  6 r  ,  6&9   and  6~e   were  then  obtained,  and  the  relations 
between  them  determined  for  the  care  where,  at  the  edge  of  the  hole 
(  r  =  a  ),  6r  ~  6re~  ?>*   The  formula  arising  from  this  analysis  is  J 

T-lien  converted  to  a  form  to  give  a  stress  concentration  factor 
referred  to  6*  ,  the  formula  becomes: 


Calculations  for  the  hole  diameters  investigated  wore  made,  a  summary 
of  which  appears  on  the  following  pajres. 
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SUMMARY  OF  STRESS  CONCENTRATION 
FACTORS  DERIVED  FROM  THE  METHOD  OF  LITRE 


3  =  0.475 


r 
& 

0° 

25° 

90° 

1.00 

5.5620 

4.7904 

1.2420 

1.25 

3.9932 

3.5633 

1.5860 

1.50 

3.3172 

3.0100 

1.5970 

1.75 

2.9691 

2.7145 

1.5437 

2.00 

2.7676 

.2.5385 

1.4850 

3.00 

2.4500 

2.2500 

1.3300 

4.00 

2.3536 

2.1585 

1.2610 

5.00 

2.3386 

2.1178 

1.2264 

10.00 

2.2582 

2.0650 

1.1778 

a 

=  0.950 

1.00 

7.2482 

6.3052 

1.9682 

1.25 

5.2497 

4.7243 

2.3077 

1.50 

4.3796 

4.0041 

2.2770 

1.75 

3.9275 

3.6163 

2.1853 

2.00 

3.6639 

3.3839 

2.0963 

3.00 

3.2445 

3.0000 

1.8757 

4.00 

3.1158 

2.8773 

1.7804 

5.00 

3.0595 

2.8226 

1.7331 

a 

=  1.250 

1.00 

8.8924 

7.7822 

2.6764 

1.25 

6.4750 

5.8564 

3.0114 

1.50 

5.4154 

4.9733 

2.9402 

1.75 

4.8620 

4.4957 

2.8110 

2.00 

4.5380 

4.2084 

2.5926 

3.00 

4.0192 

3.7314 

2.4078 

4.00 

3.8591 

3.5783 

2.2869 

5.00 

3.7886 

3.5097 

2.2272 

a 

=   1.500 

1 

1.00 

10.6052 

9.3209 

3.4U0 

1.25 

7.7514 

7.0357 

3.7444 

1.50 

6.4945 

5.9831 

3.6309 

1.75 

5.8355 

5.4117 

3.4627 

2.00 

5.4484 

5.0671 

3.3136 

3.00 

4.8264 

4.4934 

2.9620 

4.00 

4.6333 

4.3085 

2.8145 

5.00 

4.5482 

4.2250 

2.7418 
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r 

— • 

a 

1.00 
1.25 
1.50 

1.75 
2.00 
3.00 
A.  00 

5.00 


SUMMARY  OF  STRESS  CONCENTRATION 
FACTORS  DERIVED  FROM  THE  METHOD  OF  LURE 

a  5  1.7S0 


e 

0° 

25° 

90° 

12.6294 

11.1392 

4.2854 

9.2597 

8.4294 

4.6105 

7.7697 

7.1763 

4.4471 

6.9860 

6.4943 

4.2328 

6.5245 

6.0821 

4.0473 

5.7801 

5.3938 

3.6169 

5.5482 

5.1713 

3.4378 

5.4458 

5.0715 

3.3498 
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APPENDIX  C 

Stress  Concentration  Factory 

Derived  Ttojl   Flat  Plate  Theory 

Wang,  in  (7),  derives  the  stress  function  for  a  flat  plate  subjected 
to  the  action  of  a  uniform  tensile  stress  "S"  in  the  "x"  direction  as 
being;  fi~~Z$y     • 

In  terms  of  cylindrical  coordinates,  this  becomes  ^f,- 4"5rv  ~co*^  '; 
6  being  measured  from  the  "xn  axis.  From  this,  an  expression  for  cTe  is 
derived,  being*  6"0  -  \    T l   —  x"  5  ( I  ~co*  Z©)  . 

When  a  small  circular  hole  is  drilled  in  the  plate,  the  expression 
for  the  stress  function  becomes; 

%=  (C,  r^L  r  4- C2rx+Cj  lo3  r  tC4)  +(Csr2 t  Cfer  *-*?£  +  Cf)  cosZa 

When  the  boundary  conditions  are  fulfilled,  and  the  unknown  coeff- 
icients are  calculated,  the  resulting  expression  for  6&  is  derived  as 


being,  6e  =  '£:(\  +  f*)-'i(\  +  fyr)c-osZ 

By  superposing  a  tensile  stress  of  amount  2S  at  90°  away  from  d=0t 
a  condition  approaching  that  existing  in  a  thin  Bhell  under  internal 
pressure  may  be  approximated j  i.e.,  6,  ~%t  "  $    j  ^2.~  *t~*  '  %-S 

This  leads  to  an  expression  for  6e  under  the  action  of  a  uniform 
biaxial  tensile  stress  off  6e »  ^- ( If  f*J  +  -£-(/  +  ^r) c oS  Z6- 

This  equation,  when  divided  by  "5"  (the  field  stress  in  the  axial 
direction),  results  in  the  desired  expression  for  the  stress  concentration 
factor*  SCF-  f£  -  fp  =  iff  +  ■££)  +  ±( I  +  ^prjcoj  26. 

A  tabular  summary  of  values  of  SCF  for  various--  ratios  at  selected 


o- 


angles  appears  on  the  following  page. 
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SUMMARY  OF  STRESS  CONCENTRATION  FACTORS 
DERIVED  FROM  FLAT  PLATE  THEORY 

e 

r  0°  25°  90° 


a 
1.00 
1.25 
1.50 
1.75 
2.00 
3.00 
4.00 
5.00 
10.00 


5.000 

4.286 

1.000 

3.574 

3.176 

1.346 

2.963 

2.686 

1.370 

2.650 

2.4U 

1.330 

2.469 

2.257 

1.281 

2.185 

2.000 

1.148 

2.100 

1.920 

1.088 

2.062 

1.883 

1.058 

2.015 

1.836 

1.015 
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In  a  similar  manner,  it  is  possible  to  derive  an  expression 
for  ^  =  |£»  SCff   which  is; 

A  table  of  selected  values  of  this  equation  appears  below: 

e 

r 
a 

1.00 

1.25 

1.50 

1.75 
2.00 
3.00 
4.00 
5.00 
10.00 


0° 

25° 

90° 

0.000 

0.000 

0.000 

0.706 

0.646 

0.374 

0.926 

0.893 

0.741 

1.003 

1.007 

1.017 

1.031 

1.064 

1.219 

1.037 

1.142 

1.630 

1.025 

1.161 

1.787 

1.019 

1.170 

1.861 

1.005 

1.177 

1.965 
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APPENDIX  D 
SAMPLE  COMPUTATIONS  TOR  A  TYPICAL  TEST 

On  the  five  data  sheets  which  follow  (  pages  73  through  77  )■  the 
strain  indicator  readings  for  a  typical  test  -  that  for  the  1.5"  radius 
hole  -  are  shown.   Initial  readings  were  taken  at  200  psig  internal 
pressure,  and  subsequent  strain  readings  were  taken  for  each  strain 
gage  at  intervals  of  100  psi  up  to  and  including  700  psig. 

These  readings  were  in  turn  plotted  for  every  gage  and  the  linear 
portions  of  the  plots  were  used  to  determine  the  "strain  sensitivity11 
in  micro-inches  per  inch  per  psi  of  pressure.  K   tynical  exanmle  of 
this  type  of  plot  aor>ears  on  t>age  7*. 
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The  strain  sensitivities  thus  found  are  truly  indicative  ef 
the  state  of  strain  within  the  body,  and  following  eonventional 
reduction  formmlae,  can  be  reduced  to  a  stress  sensitivity  - 
psi  (stress)  per  psi  (pressure). 

The  following  reduction  formulae  from1  {11)  were  used  to 
compute  the  stress  sensitivities t 

For  biaxial  gages: 

C,  m      E   (€,  +z/££) 

1  »j/2 

&z   -  E  {-*€,  ♦  ej 

where : 

^  ■  the  algebraically  larger  principal  stress 

O^  •  the  algebraically  smaller  principal  stress 

€,   ■  the  algebraically  larger  principal  strain 

£j   -  the  algebraically  smaller  principal  strain 
E  ■  Young's  Modulus 

f   ■  Poisson's  Ratio 
For  rosettes: 

cT,  ■Ef€cu+   6C      *         l         nA€a,«  Oa+  J2^-  (^*^c)f    f 

^2(  i  -zo      Tnrr^iy *- J  J 

where : 

C,  -  the  algebraically  larger  principal  stress 

O^  ■  the  algebraically  smaller  principal  stress 
€a+  6C  -  strains  measured  by  the  perpendicular  gage  legs 

Gb    ■  strain  measured  by  the  h$°   gage  leg 
E  "  Young's  Modulus 

~&   -  Poisson's  Ratio 


,79r 


These  formulae  were  adapted  to  computation  forme,  and  as 
shown  on  the  following  pages,  were  used  to  compute  the  stress 
sensitivities  for  each  gage  location.  In  addition  to  the  con- 
ventional reduction  of  strain  rosette  data  to  principal  stress 
sensitivities,  an  additional  reduction  was  made  by  omitting  the 
U£°  leg  and  treating  the  gage  as  a  biaxial  gage,  thus  giving 
results  in  terms  of  tangential  (  Q  )   and  radial  (  r  )  stress 
sensitivities  as  related  to  the  hole0 

In  order  to  key  the  hole  instrumentation  to  the  data  reduction 
sheet,  the  following  notation  was  employed: 
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In  order  to  compare  these  stress  sensitivities  to  a  common  base, 
it  was  decided  to  use  the  theoretical  axial  stress  obtaining  in  a 
cylinder  under  pressure.  For  such  a  vessel  with  walls  of  finite  thickness, 
Seely  and  Smith  in  (12)  give  the  following  formulae  for  the  Lame  (  thick- 
wall  )  stresses: 


p  ft  > 


R/~  R] 
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6<f>: 


R*-R* 

R?  -  R? 


Equations  applicable  to  a 
cylinder  with  internally 
applied  pressure. 


By  letting  K=  %     \  »  "the  above  equations  become: 
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Then,  by  dividing  both  sides  of  the  equations  by  the  pressure 
P,  the  theoretical  stress  sensitivities  are  formed,  being  psi  (stress) 
per  psi  (pressure),  which  aret 

*--  + 

Cfe,  »  K  ♦  1 

(562.   m     K 

For  the  geometry  of  the  model  under  test  (  R^  .  7  ^  "|  j  ■  7  31"  ), 
K  ■  19c78.  Thus,  the  theoretical  axial  stress  sensitivity,  psi/psi 
equals  K/2  or  9.89o 

The  tangential  (  ©  )  and  radial  (  r  )  stresses  (  in  terms  of  a 
polar  coordinate  system  centered  on  the  hole  )  were  computed  by  aver- 
aging the  inside  and  outside  stresses  for  each  gage  location  to  ob- 
tain a  mean  stress,  which  was  then  compared  to  the  theoretical  axial 
stress.  Computations  for  this  phase  appear  on  the  following  six  pages © 
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One  of  the  general  theories  of  failure  which  seems  to  agree  well 
with  experimental  data  for  a  ductile  material  is  the  maximum  distortion 
energy  theory  attributed  to  Huber,  Hencky,  and  Von  Mises.  For  a  body 
in  a  condition  of  plane  strain,  Nadai  (13)  writes  this  as: 


♦  cr/  -  crtoi 


where<5,  and  (f.  are  principal  stresses . 
We  can  use  the  Lame  stresses  previously  determined  (pages  87  and 
88)  to  compute  a  theoretical  Hencky-Von  Mises  field  stress  sensitivity 
in  a  cylinder  for  comparison  te  the  Hencky-Von  Mises  stress  sensitiv- 
ities calculated  from  measured  data: 

Inside  Surface 


Outside  Surface 


cr* 


20,78 
9c89 
U31.808 
97.812 
529 .620 
205.  5U; 
32U.106 

18.003  psi/psi 

Hencky=Von  Mises  stress  sensitivities  were  computed  from  measured 
data  and  compared  to  the  theoretical  values  for  each  gage  location. 
Calculations  for  this  step  appear  on  pages  96  through  101. 
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19o78 

9.89 

391. 2U8 

97 « 812 
U89.060 
195. 62U 
293.U36 

17.130  psi/psi 
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One  method  of  calculating  the  stress  concentration  around  a 

reinforced  hole  is  found  in  (14-),  which  was  used  by  reference  (15) 

to  estimate  the  increased  stresses  around  apertures  in  the  hull  of 

the  "Aluminaut" .  The  following  calculations  represent  this  approach 

to  the  problem  of  determining  ;r— | 

1.  Area  of  reinforcement  (  from  reference  (H),  page  87  ) 
F-  P?  -3)(.^S  -3  7S) 
(3)  (.3  IS) 

F  =  0-23 

&  -   ■    f'  5  (  from  curve  ) 

K^  =>  2.  4-S  (  (U),  page  39  ) 

K  r  =  B  (  K^  -  /)  +  I 

Mow,  the  stress  distribution  in  terms  of  polar  coordinates  in 
an  unholed  plate  in  a  biaxial  stress  field  of  S  and  2S  may  be  taken 
as  in  (7):  . 

If  we  let  ^jX  =  S,  divide  both  sides  by  ^x,  and  multiply  by  K^  , 
we  will  have  an  expression  for  the  stress  concentration  factor  desired, 

scf  =  kk  |f  =  /.on  (j-i  cwjej 

For  the  3.0"  diameter  hole,  this  becomes: 

0  0°  30°  60°        90° 

SCF  4.03         3.51         2.51        2.04. 

Finally,  following  Murphy's  method  of  construction  of  the  Mohr's 
circle  for  strain  as  described  in  (16),  the  direction  of  principal 
strains  (  hence  stresses  )  were  determined  for  each  rosette  location 
on  the  cylindrical  surface. 


-103: 


APPENDIX  E 
DETAILED  DESCRIPTION  OF  TEST  APPARATUS 

The  test  apparatus  consisted  of  a  pressure  vessel,  pump,  flexible 
tubing,  pressure  gage,  and  equipment  required  for  the  taking  of  strain 
gage  data.  The  details  of  the  pressure  vessel  are  shown  in  Figures  11 
through  H,  pages  106  through  109. 

The  pressure  vessel  was  constructed  of  HY-80  steel  in  accordance 
with  a  basic  design  drawn  up  by  Mr.  John  Pulos  of  the  David  Taylor 
Model  Basin,  The  fabrication  sequence  was  as  follows: 

1)  A  1"  thick  HY-80  plate  was  rolled  to  a  14."  I.D.  and  closed 
with  a  longitudinal  weld  to  form  a  cylinder. 

2)  The  cylinder  was  stress  relieved  at  1100°F  for  two  hours. 

3)  The  cylinder  was  machined  to  the  final  I.D.  and  O.D. 
4.)  The  radially  oriented  hole  was  machined  out  and  the 

reinforcement  plug  welded  in. 

5)  The  cylinder  was  stress  relieved  at  1100°F  for  two  hours. 

6)  The  pressure  vessel  end  rings  were  welded  on. 

The  pressure  vessel  was  instrumented  with  electrical  resistance 
foil  strain  gages  as  shown  in  Figure  15.  In  addition,  gages  were 
placed  on  the  inside  periphery  of  the  hole  in  the  reinforcing  plug. 

The  biaxial  strain  gages  located  on  the  0°  and  90°  legs  were  made 

up  of  two  Budd  C6-111  foil  gages  placed  perpendicularly  to  each  other. 

The  rosette  strain  gages  on  the  25°  leg  were  Baldwin  SR-4.  type  FABR 
12-12  foil  gages  (for  gage  locations,  see  Figure  A,  page  7). 
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The  hole  periphery  strain  gages  for  hole  diameters  of  .95",  1.90" 

and  2.50"  were  Budd  C6-111  foil  raras  Th.   k-i     ,   ■. 

°u  Sages.  Ihe  hole  periphery  strain 

gages  for  hole  diameters  of  3.00"  and  3.50"  were  Baldwin  SR-4 

•  type  FABR  12-12  foil  gapes. 

All  strain  gages  on  the  0°,  25°  and  90"  legs  were  applied  with 
HYS0L  epoxy  cement,  following  the  manufacturer's  instructions.  All 
strain  gages  except  those  on  the  hole  periphery  were  covered  with 
epoxy  cement  for  their  protection  after  being  wired  to  the  strain 
indicator  and  tested  for  efficiency  of  bonding  by  pressing  with  a 
pencil  eraser  while  watching  for  deflections  of  the  strain  indicator 
needle.  Strain  gages  on  the  hole  periphery  were  applied  with  Eastman 
910  adhesive  and  left  uncovered  since  they  were  only  used  for  a 
short  period  of  time  and  had  to  be  replaced  for  each  hole  size. 

The  dummy  reference  strain  gages  were  attached  to  a  block  of 
HY-80  steel  which  was  located  inside  the  pressure  vessel  to  keep 
all  strain  gages  at  or  near  the  same  temperature  while  strain 
readings  were  being  taken. 

The  actual  arrangement  of  the  test  apparatus  is  illustrated 
in  Figures  2  and  3,  pages  5  and  6. 

A  detailed  list  of  apparatus  appears  on  pare  105. 
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LIST  OF  EQUIPMENT  ASSOCIATED  WITH  PROJECT 

A.  STRAIN  INDICATOR: 

Baldwin  -  Lima  -  Hamilton,  Type  N 
B#  STRAIN  GACE  SWITCH  BOX: 

Type  186-C,  4&  position 

David  W.  Taylor  Model  Basin 

Washington,  D.C. 

C.  HYDRAULIC  PUMP: 

Blackhawk  Hydraulic  "Porto-Power"  Jack 
Type  P  -  76,  0  -  20,000  psi,  hand  operated 

D.  FLEXIBLE  TUBING: 

Blackhawk  "Porto-Power"  flexible  hose,  Z  -  864 

SAE  100R1,  3/16"  wire  reinforced 

Equipped  with  ZH  -  630  Bantam  "SPEE  -  D  -  COUPLER" 

E.  PRESSURE  GAGE: 

Ashcroft  bourdon  tube  type 

0  -  1,000  psi  in  10  psi  subdivisions 

8-p'  dial  face 

F.  HYDRAULIC  OIL: 

Lubricating  oil,  general  purpose 
Navy  symbol  304.2;  MIL  -  L  -  15016A 
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RIAL:      HY     RO   STEEL 
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PRESSURE     VESSEL      CLOSURE     PLUG 
FIGURE     |Z 
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MATERIAL:  MEDIUM     STEEL 
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TAP  i-     -7,^  NF,   I"  DEE  P 
WHOLES    EQUALLY    oPACED 
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ON     I5.5*PIA.   &.C     TO     MATCH 
HOLE  S    IN    END     RinjG 


PRESSURE     VESSEL      COVER     PLATE. 


FIGURE.     14 
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FIGURE  15 
Photograph  of  Pressure  Vessel  Instrumentation 


-110- 


FIGURE  16 
Photograph  of  Wiring  Set-Up 
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APPENDIX  F 
MISCELLANEOUS  PROCEDURES  AND 

VARIATIONS  IN  BASIC  TECHNIQUES 


SECTION  I 
INFORMATION  REGARDING  THE  ATTACHMENT 
OF  METAL  FOIL  STRAIN  GAGES 


SECTION  II 

REPLACEMENT  OF  CLOSURE  PLUG  O-RING 

WITH  A  LOW  ADHESION  SEALING  COMPOUND 


SECTION  III 
REMOVAL  OF  REINFORCEMENT 


SECTION  IV 
DETECTION  OF  A  PRESSURE  SENSITIVE  GAGE 
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SECTION  I 
INFORMATION  REGARDING  THE  ATTACHMENT 
OF  METAL  FOIL  STRAIN  GAGES 

Two  different  cement  kits  were  used  in  the  attachment  of  the 
bakelite-baoked  metal  foil  strain  gages  (Budd  C6111  and  Baldwin  SR4 
FABR  1^-12  types)  used  to  instrument  the  periphery  of  the  hole.  Both 
of  these  kits  are  based  on  the  use  of  Eastman  910  adhesive,  a  cyan©- 
acrylate  monomer  which  can  transform  almost  instantaneously  from  a 
free-flowing  liquid  to  a  rigid  plastic,  forming  strong  bonds  with 
almost  any  material. 

Of  the  two  kits  used,  the  one  purchased  from  the  Baldwin-Lima- 
Hamilton  Corporation  gave  unsatisfactory  results,  although  the  instruc- 
tions contained  therein  were  followed  explicitly.  On  the  other  hand, 
the  GA-1  contact  cement  kit  supplied  by  the  Budd-Tatnall  Company  gave 
extremely  good  results.  The  procedure  set  forth  in  the  Budd  Company 
instructions  was  modified,  however,  on  the  basis  of  information  received 
from  the  David  Taylor  Model  Basin  Personnel,  For  possible  future  refer- 
ence, the  method  which  the  authors  found  to  be  satisfactory  is  outlined 
below  s 

A,  Surface  Preparation 

1,  Remove  burrs,  tool  marks,  etc.,  from  the  surface,  A  small 
grinding  wheel  followed  by  No,  2  emery  paper  seems  suitable, 

2,  Using  cotton  applicators,  clean  the  area  thoroughly  with 
benzene  until  the  swab  tips  are  completely  clean.  Follow  up 
with  acetone  cleaning  in  the  same  manner. 


-114- 


3.  Apply  Budd  0A-1B  neutralizer  with  a  swab,  let  stand  a 
few  seconds,  then  remove  with  dry  swabs  (GA-1B  is  a  solution 
of  ammonia  in  water;  it  has  been  found  that  Eastman  910  sets 
up  best  when  the  environment  is  very  slightly  basic). 
£•  Apply  GA-1A  accelerator  to  the  surface.  This  must  be 
allowed  to  dry  for  at  least  four  minutes  before  placing  the 
gage  in  position. 

B.  Gage  Preparation 

1.  Apply  cellophane  tape  ("Scotch  Tape")  over  surface  of 
the  gage  and  peel  off  any  temporary  backing.  Clean  gage 
back  with  acetone. 

C.  Installation 

1.  Apply  a  drop  of  910  cement  to  back  of  gage.  Spread  out 
evenly,  using  the  stick  of  a  cotton  swab,  but  not  exerting 
any  direct  pressure  on  the  gage. 

2.  Apply  gage  to  the  surface,  and  press  firmly  into  position 
with  finger.  Avoid  getting  cement  on  the  fingers,  as  they 
will  adhere  firmly  to  the  surface.  After  one  minute,  the 
carrier  tape  may  be  peeled  slowly  off  of  the  gage,  and  wire 
connections  may  be  made. 
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SECTION  II 

REPLACEMENT  OF  CLOSURE  PLUG  O-RING 

WITH  A  LOW  ADHESION  SEALING  COMPOUND 

For  the  last  two  tests  in  the  series,  the  C-ring  type  of  seal  on 
the  closure  plug  would  not  have  been  satisfactory  because  of  the  geo- 
metric limitations  imposed  on  the  plug  and  the  lack  of  sufficient  fay- 
ing surface  against  which  an  O-ring  must  bear.  Following  the  advice  of 
David  Taylor  Model  Basin  personnel,  a  low  adhesion  sealing  compound  was 
employed  with  excellent  results.  This  compound  is  made  by  the  Products 
Research  Company  of  Gloucester  City,  New  Jersey,  and  is  designated  as 
Sealant  PR-1321. 

The  sealant  is  a  synthetic  rubber-base,  Thiokol  liquid  polymer, 
compounded  into  a  highly  thixotropic  red  paste  which  may  be  applied  to 
the  surfaces  to  be  sealed  by  a  spatula.  It  cures  to  a  solid  rubber 
which  seals  to  the  extent  that  no  leakage  at  all  was  observed  at  any 
time  during  the  tests. 

The  procedure  which  the  authors  found  quite  satisfactory  is  as 
follows : 

a.  Wash  the  faying  surfaces  thoroughly  with  benzine,  followed 
by  acetone, 

b.  Thoroughly  mix  the  sealing  compound;  10  parts  by  weight  of 
base  compound  to  1  of  accelerator, 

c.  Apply  sealant  to  one  of  the  faying  surfaces,  then  join  them. 
Apply  moderate  to  heavy  pressure  to  cause  the  compound  to 
flow  into  the  surface  irregularities.  Let  stand  under  this 
pressure  for  24.  hours  for  curing. 
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For  the  tests  in  which  this  compound  was  used,  the  sealant  was 
forced  into  the  surfaces  by  tightening  up  the  closure  plug  bolt  and 
strongback  arrangement.  During  the  actual  testing  period,  the  closure 
plug  nut  could  be  backed  off  with  no  danger  of  leakage.  When  subsequent 
disassembly  of  the  model  was  performed,  it  was  found  necessary  to 
strike  the  plug  smartly  with  a  mallet  in  order  to  break  the  seal  thus 
formed. 
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SECTION  III 
REKOVAL  OF  REINFORCEMENT 

Successive  enlargements  of  the  hole  (  or  removal  of  reinforcement  ) 
proved  to  be  a  more  tedious  and  difficult  task  than  had  been  originally 
conceived.  A  cradle  of  some  sort  was  needed  in  which  to  mount  the  model 
during  machining,  so  an  old  engine  assembly  bed  was  modified  to  do  the 
job.  This  permitted  the  model  to  be  clamped  securely  to  the  bed  and  the 
bed  in  turn  bolted  to  the  base  of  the  Wiener  drill  press  in  the  machine 
shop. 

The  Wiener  press  which  was  used  to  machine  away  the  reinforcement 
cannot  be  considered  adequate  for  a  Job  of  this  type.  When  the  machine 
was  being  used  to  bore  out  the  hole  in  the  model,  the  entire  "fixed" 
head  of  the  drive  and  spindle  could  be  observed  describing  an  elliptical 
path  180°  out  of  phase  with  the  cutting  tool.  Only  by  taking  small  cuts 
at  low  feed  rates  could  reasonably  circular  holes  be  machined. 

New  tools  were  purchased  by  the  Institute  to  assist  in  the  job  of 

t 

reinforcement  removal;  a  Chandler-Duplex  Model  J  Combined  Boring  and 
Facing  Tool  Head,  and  a  "J"  set  of  Bokum  Boring  Bars.  The  Chandler-Duplex 
"J"  is  a  versatile  tool  head  which  permits  boring  out  a  hole  up  to  6W  in 
diameter,  or  facing  a  like  area.  (  The  latter  feature  was  not  employed 
by  the  authors.)  The  set-up  used  can  be  seen  in  Figure  18.  Because  of  the 
lack  of  flexural  rigidity  of  the  machine  in  which  it  was  used  as  well  as 
the  toughness  of  the  HY-80  steel  which  was  being  machined,  the  maximum 
cut  that  could  be  taken  was  .050"  on  a  diameter.  (  Even  this  modest  cut 
produced  occasional  chattering. )  It  is  felt  that  the  Model  J  head  did  not 
receive  a  fair  trial  in  this  respect,  and  that  when  used  on  softer  material 
or  in  a  better  machine,  it  will  perform  very  creditably. 
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FIGURE  18 
Photograph  of  Reinforcement  Removal 
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SECTION  IV 
DETECTION  OF  A  PRESSURE  SENSITIVE  GAGE 

While  the  tests  were  being  performed  with  concurrent  reduction  of 
strain  data,  it  was  found  that  all  of  the  strain  gage  rosettes  on  the 
inside  25°  leg  were  producing  readings  far  in  excess  of  what  should 
normally  have  been  expected.  As  the  model  was  filled  with  oil  for  testing, 
access  to  the  suspect  "dummy"  gage  associated  with  this  rosette  group 
was  impractical.  In  order  to  verify  the  fact  that  this  dummy  was  pressure 
sensitive,  it  was  connected  up  as  an  "active"  gage  to  the  strain  indica- 
tor while  a  new  rosette  was  mounted  on  a  solid  metal  block  and  connect- 
ed up  as  the  dummy.  When  pressure  was  applied  to  the  model,  deflection 
of  the  indicator  needle  immediately  revealed  that  the  dummy  gage  in 
question  was  indeed  pressure  sensitive. 

Using  the  newly  mounted  exterior  dummy,  the  other  rosette  dummy 
gages  were  checked  and  found  to  be  satisfactory  (  no  changes  in  strain 
for  changes  in  pressure  were  observable  )#  The  test  dummy  was  removed 
from  the  circuit,  and  a  regular  series  of  strain  readings  were  taken 
on  the  inside  rosettes  using  first  the  pressure-sensitive  dummy  and 
then  a  good  dummy  in  the  circuit.  Strain  sensitivity  plots  were  made 
and  both  found  to  be  linear,  so  a  correction  factor  was  computed  for 
each  gage  element  and  applied  to  the  readings  of  the  affected  gages 
for  the  previous  test  (  1.9"  diameter  hole  ). 

Upon  later  disassembly  of  the  model  for  machining,  the  block  of 
metal  with  the  dummy  gages  mounted  theron  was  inspected.  It  was  found 
that  the  suspected  dummy  rosette  had  indeed  become  unbonded,  thus 
causing  erroneous  readings  to  result.  The  lack  of  bond  was  quite 
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apparent  for,  when  poked  with  a  rubber  pencil  eraser,  violent 
deflections  of  the  strain  indicator  needle  resulted. 

It  is  felt  that  the  practice  of  placing  all  dummy  gages  in  the 
interior  of  the  model  is  not  entirely  satisfactory.  Although  it 
prevents  inadvertent  damage  to  them  and  insures  that  they  will  not 
be  subject  to  temperature  fluctuations,  it  makes  access  to  them 
rather  difficult. 

It  is  believed  that  a  better  arrangement  would  be  to  have  the 
dummies  associated  with  the  interior  gages  mounted  on  a  gage  block 
inside  the  model,  and  those  associated  with  outside  instrumentation 
left  exterior  to  the  model.  In  this  way,  should  any  bad  readings 
occur  on  at  least  the  exterior  gages,  the  simple  and  expedient 
"poke  test"  with  a  pencil  tip  eraser  could  be  immediately  used. 
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